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Molecular motors and cytoskeletal filaments work collectively most of the time under opposing 
forces. This opposing force may be due to cargo carried by motors or resistance coming from the 
cell membrane pressing against the cytoskeletal filaments. Some recent studies have shown that 
the collective maximum force (stall force) generated by multiple cytoskeletal filaments or molecular 
motors may not always be just a simple sum of the stall forces of the individual filaments or motors. 

To understand this excess or deficit in the collective force, we study a broad class of models of both 
cytoskeletal filaments and molecular motors. We argue that the stall force generated by a group of 
filaments or motors is additive, that is, the stall force of N number of filaments (motors) is N times 
the stall force of one filament (motor), when the system is in equilibrium at stall. Conversely, we 
show that this additive property typically does not hold true when the system is not at equilibrium at 
stall. We thus present a novel and unified understanding of the existing models exhibiting such non- 
addivity, and generalise our arguments by developing new models that demonstrate this phenomena. 

We also propose a quantity similar to thermodynamic efficiency to easily predict this deviation from 
stall-force additivity for filament and motor collectives. 


I. INTRODUCTION 

Molecular motors (such as kinesin, dynein and myosin) 
and cytoskeletal filaments (such as actin and micro¬ 
tubule) are abundantly present in eukaryotic cells and are 
responsible for important cellular functions. Cytoskele¬ 
tal hlaments give structural stability to the cell and act 
as tracks along which molecular motors move and facili¬ 
tate intra-cellular transportation M- Many researchers 
have studied the dynamics of single hlaments and sin¬ 
gle motors in great detail using experimental and the¬ 
oretical tools [TH^. However, inside a cell, molecular 
motors or cytoskeletal hlaments work collectively most 
of the time to perform their biological tasks. For ex¬ 
ample, the dynamics of multiple actin hlaments is re¬ 
sponsible for lamellipodial protrusions during crawling 
of cells [U . Similarly, microtubules work collectively 
to bring about segregation of chromosomes during mito¬ 
sis [1^. Many dynein motors attach to molecular cargo 
and generate force for cellular transportation [ 23 , 1^ . 
while myosin motors primarily work together to generate 
forces in stress hbres and muscle tissues. Hence, study 
of such systems, which are involved in a wide range of 
biological processes, requires understanding of the gener¬ 
ation of collective force by multiple cytoskeletal hlaments 
or motors [23 - [^ . 

To understand collective force generation by polymer¬ 
izing biohlaments, researchers typically resort to in vitro 
experiments in which biohlaments are polymerised either 
in the form of bundles or branched sheets against a bar¬ 
rier, which resists their motion by producing reaction 
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forces [H, l39l - l43| . One important focus of these stud¬ 
ies is a quantity called the stall force, which is dehned as 
the resisting force at which the mean growth velocity of 
the collection of hlaments is zero, and is the maximum 
pushing force generated by these hlaments. 

A number of experiments on force generation by an 
assembly of biohlaments have been reported in the liter¬ 
ature [Sfi, m, Idll - li^l . In view of the overall content and 
objective of our paper we will focus only on the exper¬ 
iments related to single and parallel bundles of growing 
biohlaments. Experiments on single actin-hlament are 
relatively rare due to experimental limitations created 
by the easy buckling tendency of a single actin hlament. 
To the best of our knowledge, only one set of experiments 
has reported stall force for a single growing actin hlament 
- a value of ~ IpN [4l|. In another set of experiments 
on a few actin hlaments growing in parallel. Footer et 
al. [i^ reported a collective stall force of « I pN, which 
is very similar to the load required to stall a single hla¬ 
ment. In this experiment, the interpretation of hlament 
cooperativity while pushing together against the barrier 
is, however, complicated by the fact that individual actin 
hlaments may buckle before reaching their stall force, and 
hence, collectively the bundle may be unable to reach its 
full potential for force generation. On the other hand, 
in similar experiments performed on microtubules, re¬ 
searchers have found that the growth velocity of the hla¬ 
ment growing against an immobile barrier decreases with 
increased resistance - from 1.2/rm per minute at zero 
load to 0.2pm per minute at 4pN toforce, which is 
the putative stall force in this case |39|. In another ex¬ 
periment based on microtubules, using optical tweezers, 
Laan et al. report that stall forces of 2.7pN, 5.5pN 
and S.lpN are produced by groups of hlaments. They in- 
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terpreted this as containing one, two and three filaments, 
respectively. Hence, they concluded that stall force scales 
linearly with the number of filaments though there could 
be possible ambiguity in directly counting the number of 
filaments in the group. 

On the theoretical and computational front, there are 
a number of very detailed models for understanding the 
force-velocity dynamics of a single biofilament jSOl - l^ , as 
well as a few that try to model the dynamics of multiple 
biofilaments [1^ . Some of these theoretical 

studies have demonstrated that the stall force of multiple, 
non-interacting filaments without ATP/GTP dynamics, 
scales linearly with the number of filaments [1^ [2^ 
[13, 115 . In contrast, a few other recent papers quite 
clearly report that inclusion of ATP/GTP hydrolysis can 
lead to either enhanced or reduced cooperativity in the 
maximum force generated by multiple growing filaments; 
the stall forces need not always scale with the number of 
filaments (111. Is^-is^. In other words, the stall forces of 
individual filaments are non-additive in general, that is, 
the collective stall force produced by N number of fila¬ 
ments (denoted by /i^^) is not just a simple sum of indi¬ 
vidual stall forces of single filaments (i.e., ^ Nfs^'^). 

It is quite intriguing as to why tweaking an internal vari¬ 
able (ATP/GTP) for an individual filament without ac¬ 
tually changing any external interaction between the fil¬ 
aments can lead to this drastic change in their collective 
generation of force. 

In studies similar to those on the biofilaments, the stall 
force for motors is defined by the resisting force at which 
the average velocity of the motors is zero. The exper¬ 
iments on molecular motor force generation mostly in¬ 
volve micro-sized dielectric particles and optical tweez¬ 
ers, by which a resisting force is applied on the moving 
motors, in order to measure their velocity response as a 
function of the resisting force We briefly describe 

in the following text, a few examples that are relevant for 
our current study. Single molecule study of kinesin shows 
that kinesin is a strong molecular motor and generates 
maximum force of magnitude « 5-8pN [H, , whereas 

the stall force of certain variants of dynein is measured 
to be a relatively lower value of Ri IpN [b^l. Dyneins in 
a group are, however, reported to generate force collec¬ 
tively, something that is missing in kinesin, mainly be¬ 
cause of its lack of processivity under larger forces . A 
single headed variant of the kinesin family, KIFIA, which 
migrates along the microtubule in alternating states of 
strong attachment and incomplete detachment, produces 
a very small stall force (« 0.1 pN). However, a very recent 
experiment of tube-pulling assay on KIFIA motors has 
demonstrated extremely strong cooperative force gener¬ 
ation in these motors - ss 10-15 single headed KIFIA 
motors could indeed pull out tubes from giant unilamel¬ 
lar vesicles, which require a force around two orders of 
magnitude larger than the arithmetic sum of the individ¬ 
ual stall forces 0. 

There are a few broad classes of models present in 
the literature for understanding the force-velocity rela¬ 


tion of both single molecular moto rs fl3.1^ . as well 
as for a group of molecular motors [74l - l7^ . The multiple 
molecular motor models describe a variety of different 
biophysical scenarios such as motors elastically coupled 
to each other, motors elastically coupled to the cargo, tug 
of war between motors walking in opposite directions and 
self-exclusion interactions between motors pulling on a 
membrane tether, for processive as well as non-processive 
motors (3 ^ Specifically, Campas et al. [11] have 
shown theoretically that the collective stall force for mul¬ 
tiple motors are non-additive (/i^^ ^ Nfs^'^) in the pres¬ 
ence of attractive or repulsive interactions and can be 
manyfold larger. However, in the absence of such in¬ 
teractions, the stall forces in this model are simply ad¬ 
ditive. Also, in a series of recent papers, Casademunt 
and co-workers [H, [H, [tII using a ‘two-state’ model [ill 
for multiple interacting motors, have demonstrated that 
the motors can produce highly enhanced cooperativity in 
stall force generation, and, in particular, demonstrated 
this phenomenon for KIFIA motor 0. 

From the above discussion, it is quite clear that the 
collective force generation by motors and biofilaments 
can indeed exhibit or has the potential to exhibit highly 
diverse behaviour. As noted earlier, some studies have 
shown enhanced/reduced cooperativity in the collective 
stall force generation, while other studies hint towards 
additivity of stall forces. Hence, it would be both inter¬ 
esting and important to understand the conditions un¬ 
der which stall forces become simply additive, and con¬ 
sequently, get better insight into the circumstances un¬ 
der which the simple additivity is lost. In this paper, we 
develop a general theoretical framework to understand 
how enhanced/reduced cooperativity in collective force 
generation can appear in systems of multiple cytoskele- 
tal filaments or motors. We investigate this question by 
studying various models for these systems. From our 
case studies we conclude, quite generally, that the vio¬ 
lation of stall force additivity stems from the violation 
of the condition of detailed balance, that is, departure 
from thermodynamic equilibrium. On the other hand, 
stall forces are additive when the system of filaments or 
motors is in equilibrium at stall. The main contribution 
of this paper is to recognize the hitherto invisible thread 
of thermodynamic equilibrium linking the phenomenon 
of stall force additivity across a variety of models. 


II. COLLECTIVE STALL FORCE FOR 
MULTIPLE CYTOSKELETAL FILAMENTS AND 
MOTORS: STALL FORCES ARE ADDITIVE FOR 
EQUILIBRIUM DYNAMICS 

Inspired by the growth of cytoskeletal filaments against 
cell membranes, theorists have studied various models 
of filament dynamics against a constant applied load. 
HE HE mi He he mj. similar to the filaments, motors 
also work against load and have been modelled in several 
studies (sEI^ItE I^. Many of these filament and motor 
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models are mathematically very similar and can be de¬ 
scribed by a common framework of biased random walk 
(see FigHKa) and FiglTJb)). We consider a collection of 
rigid filaments nucleating from a fixed wall at left, while 
a resisting force / is applied upon them via a moving wall 
at right. Each filament grows and shrinks by stochastic 
addition and removal of subunits of size d (see Fig. [Ija)). 
Similarly, each motor takes a single forward or backward 
step on a fixed one-dimensional track (such as a track of 
microtubule or actin) of lattice constant d (see Fig. [Hb)). 
However, there is a key difference between the systems of 
filaments and motors: the motors cannot overtake each 
other, maintaining their sequence on the lattice, and con¬ 
sequently the leading motor alone bears the force / (Fig. 
[T{b)). The motors further obey ‘mutual exclusion’, that 
is, each lattice site can be occupied by one motor only 
when the site is empty. On the contrary, the filament- 
tips do overtake each other since the filaments grow par- 
allelly, and, therefore, any of the filaments can experience 
the applied force / (Fig. [Ha)). In these models, we mea¬ 
sure the forces in the units of ksT/d, where ks is the 
Boltzmann constant, T is the temperature, and d is the 
subunit size (or lattice size). We take ksT/d = 1 with¬ 
out losing generality. In the absence of any force, we de¬ 
note the polymerization rate (forward-hopping rate) and 
depolymerization rate (backward-hopping rate) for fila¬ 
ments (motors) as u and w respectively. In the presence 
of a resisting force /, the polymerization rates (forward¬ 
hopping rate) decreases, and the depolymerization rate 
(backward-hopping rate) increases according to the fol¬ 
lowing rules: u{f) = ue~^^ and w(/) = . Here, 

the parameter S ^0, 1] is commonly known as force dis¬ 
tribution factor [cite for motors also]. 

In the context of the models mentioned earlier, we can 
argue that the system is in equilibrium at stall. By def¬ 
inition, at the stall force, the mean velocity of the sys¬ 
tem is zero which implies that the overall monomer flux, 
in and out of the filament assembly, is zero. Since the 
monomer flux is zero, it is logical to believe that the sys¬ 
tem is in thermodynamic equilibrium at stall. Strictly 
speaking, the system is in equilibrium at stall, only if it 
is bounded in length. Unbounded systems will be freely 
diffusive at stall, which is clearly a non-equilibrium phe¬ 
nomenon. In fact, our class of systems (filament or motor 
collectives) have finite sizes for all practical purposes. In 
the case of motors, they can be thought to diffuse on 
a tilted energy landscape (see FiglHd)), which results 
in the biased random walk. This bias can be created 
by chemical potential that is linked with the ATP hy¬ 
drolysis [13, which is irreversible and hence inherently 
a non-equilibrium process. However, in the context of 
these models, the role of the ATP chemical potential is 
just to create a tilt in the energy landscape. Thus, as far 
as the translational movement of the motors against an 
applied force is concerned, we can conceptually think of 
the system to be in thermodynamic equilibrium at stall 
with respect to the translational degree of freedom - the 
chemical coordinate is simply orthogonal to the transla¬ 


tion coordinate. This seems quite analogous to the case 
of non-interacting active Brownian particles, which exert 
pressure on the confining walls similar to an ideal gas 
in equilibrium, but with a renormalised temperature due 
to the free energy consuming activity [111. Nevertheless, 
these arguments cannot be claimed to be true for every 
model for motors or filaments. In fact, in the subsequent 
sections, we show that for most models, the arguments 
advanced here break down and the systems are not in 
equilibrium at stall in general, except for certain choices 
of the parameters. 



FIG. 1. Biased random walk model for filaments and mo¬ 
tors: (a) Multiple {N > 1) filaments against a wall (see [2^ 
for a detailed study of this model). Polymerization and de¬ 
polymerization rates are u and w, when the filaments are away 
from wall (hence force-free). Note that, when more than one 
filament touches the wall simultaneously, the depolymeriza¬ 
tion rate becomes independent of the force)w), as a single 
depolymerization event cannot cause any wall movement, (b) 
Multiple motors moving on a track. Note that only the lead¬ 
ing motor bears the force /, while the trailing motors are 
experiencing force through hard core interaction, (c) Two 
filaments at the collective stall force /j ' undergoing simple 
processes of polymerization and depolymerization, (d) Two 
motors diffusing on a tilted free-energy landscape. The force 
is acting only on the leading motor, and the energy land¬ 
scape of the leading motor changes with the application of 
force. Due to hard-core interaction, a larger force(> /i^^) is 
required to stall the system of two motors. It is like moving 
uphill, with the assistance of a trailing motor, which rectifies 
the average backward motion of the leader [^ . 


Based on these arguments, if one assumes that the sys¬ 
tem is in equilibrium at stall, the tools of equilibrium 
statistical mechanics can be applied to calculate the stall 
forces. Here we show that the stall forces are additive 
for multiple filaments in the simplest model, where only 
polymerization and depolymerization processes are in¬ 
volved (see FiglHc)). We first consider the dynamics of 
a single filament under the stall force / = applied 
via the right wall. Let the wall position be x in terms of 
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the subunit size (see FiglUc)). Using equilibrium statis¬ 
tical mechanics we can write the probability distribution 
of the wall-position P{x) as below 


1 


P{x) = -e 


— /3ex _ J_ -/3(/i^^-£)x 




( 1 ) 


sition is 


P{x) = 


xi=0 


e 2 €)x ^ large 


(3) 


where Z is the partition function, and e = ln(u/i/;) is 
the free energy gained per subunit through polymeriza¬ 
tion. Note that the term e~^^‘ appears as we have a 
Gibbs ensemble in 1-dimension with fixed external com¬ 
pressive force However, this distribution P{x) 

should be independent of x, as the stall condition can 
be reached at any position x. As per Eq. [T1 this implies 
that = e = \ti{u/w). This same expression has been 
obtained in earlier theories from a detailed kinetic cal¬ 
culation [ 2 ^. Next, we consider a two-filament system 
subjected to their stall force as shown in FiglT](c). 
Let the tip-position of the trailing filament be xi, which 
is between 0 and the wall position x. The probability 
distribution of the wall-position, if the system is in ther¬ 
modynamic equilibrium, is: 


P{x) = ^ 

0 

^ large X 


( 2 ) 


A factor of 2 appears on the right hand side, since there 
could be two equally likely situations - either the top- 
filament or the bottom-filament can be the leader (see 
Fig- me)). Again, P{x) is expected to be independent of 
X, implying that = 2e = 2/i^^ This argument can be 
easily extended to > 2, and thus for this 

simple model. Some arguments, based on detailed bal¬ 
ance criterion, have also been given in Refs. [ 2 ^ and [ 2 ^ 
to show similar result of oc N for their respective 
models on cytoskeletal filaments. However, we demon¬ 
strate that a simple calculation based on elementary sta¬ 
tistical mechanics, for essentially kinetic processes, leads 
to the same conclusion. 

We now develop similar arguments to show the additiv¬ 
ity of stall forces for multiple motors. The forward and 
backward hopping processes for motors can be viewed 
as random walks on a tilted free-energy landscape (see 
FigHKd)). In this case, x and xi should be interpreted 
as the positions of the leading and the trailing motors 
respectively, for a two-motor system. The free energy 
‘released’ per unit step by going downhill on the free- 
energy landscape is e (equivalent to the polymerization 
energy). For a single motor under stall force, we can 
write the same equation as before (Eq. [T]) for the proba¬ 
bility distribution of the leading motor’s position, P{x), 
by recognizing the system to be in equilibrium at stall. 
However, for two motors at stall, the probability distri¬ 
bution of the leader’s position is somewhat different from 
the previous case of the filaments, since motors cannot 
overtake each other. The distribution of the leader’s po- 


Again, using the argument that P{x) should be indepen- 

( 2 ) 

dent of X, we get back the force additivity : fs=2e = 

2/i'^ 

In summary, we show that the stall forces are addi¬ 
tive for the simple models considered here. This demon¬ 
stration hinges on the recognition that the systems are 
in equilibrium at stall, which may not be true in many 
biological situations. In fact, we show in the next sec¬ 
tions that there are many classes of models for which 
the equilibrium description is simply not feasible, and 
consequently more complex models have interesting im¬ 
plications. 


III. STALL FORCES ARE NON-ADDITIVE FOR 
BIOLOGICALLY RELEVANT 
NON-EQUILIBRIUM MODELS 

In this section, we present several case studies to show 
that the force inequality (/i^^ ^ Nfs^'^) is true in general 
for stall dynamics departing from equilibrium - however, 
for certain combinations of kinetic rates the relationship 
= Nfs^'^ can be retrieved. We begin by analyzing 
various models of cytoskeletal filaments. 


A. A random hydrolysis model for cytoskeletal 
filaments 

In cytoskeletal filaments (such as, microtubules and 
actin filaments), subunits are typically bound to 
ATP/GTP molecules. When the subunits are connected 
to the filaments, the ATP/GTP molecules release phos¬ 
phate and convert to ADP/GDP in a process known as 
ATP/GTP hydrolysis [1, The ADP/GDP-bound 
monomers have much higher depolymerization rates com¬ 
pared to ATP/GTP-bond monomers Due to this 

heterogeneity the cytoskeletal filaments exhibit interest¬ 
ing properties such as ‘dynamic instability’ [s^- The 
dynamics of the cytoskeletal filaments have been theo¬ 
retically studied by many researchers using the “random 
hydrolysis” model Here, we focus on the simpli- 

fied model of this process as discussed in Refs. [^ . 
where we neglect complexities related to biofilaments 
such as structural details, mechanical flexibility, possi¬ 
ble interplay between mechanical forces and hydrolysis 
events and the possible multi-stage nature of the hydrol¬ 
ysis process [13, HI, [H, |^, [111. As shown in Fig. [2Ka) , 
we consider multiple filaments undergoing random hy¬ 
drolysis and growing against a wall held by a constant 
opposing force /. In the model, each monomer can be in 
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FIG. 2. Random hydrolysis model with three filaments, 
where individual monomer switches from T to D unidirection- 
ally and randomly. Different processes are shown by arrows. 
ko = 3.2 c = 100/rM, wt = 24s“^, and r = 0.2s“^. 

(a) (red curve) and a (blue curve) versus wd- Note that 
both a and are zero only at wt = wd and they are 

highly correlated in sign, (b) The fluxes per filament (defined 
in the text), at stall, for T and D monomers as a function 
of the filament number N for wd = 290s“^ (a > 0). (c) 
The collective stall force per filament /N) against N, 
the number of filaments. For a > 0, wd = 290s“^ and for 
a < 0 , Wd = 5s“^ 


two states: T (ATP/GTP-bound), and D (ADP/GDP- 
bound). Only T monomers bind to the filaments with a 
rate u{f) = ue^ (next to the wall) or u (away from the 
wall). The rate u is proportional to the concentration 
(c) of T monomers and is defined as tt = k^c. The de¬ 
polymerization occurs with a rate wt if the tip-monomer 
is T, and wd if it is D. For simplicity we assume that 
there is no force dependence on the off-rates wt and wd 
(i.e., (5 = 1). Hydrolysis (T to D conversion) happens 
randomly on any T subunit in space with a rate r. Note 
that the conversion T—:>D is irreversible, as it is not bal¬ 
anced by a reverse conversion, which makes the dynamics 
inherently non-equilibrium. The exact analytical result 
for the stall force is not known for such a detailed 
model. Instead we numerically simulate the model us- 

see Appendix A) with 


for microtubules 


ing the Gillespie algorithm 1 
experimentally known rates j 
and actin filaments. 

Before proceeding to discuss the issue of stall force ad¬ 
ditivity, we define a parameter, the “force deviation”: 


^(N) ^ fiN) _ ^ jil) 


(4) 


This parameter represents the excess/deficit of the collec¬ 
tive stall force generated by N filaments, as compared to 
N times the force generated by a single filament. So the 
deviation 0 implies the violation of force equality. 

As noted before, the random hydrolysis model is a 


non-equilibrium model by construction. To understand 
the implications of this non-equilibrium nature, we first 
look at the energetics associated with the polyineriza- 
tion/depolymerization processes. A growing filament 
clearly performs work against an applied load / through 
polymerization. The work done for addition of one sub¬ 
unit to the filament is simply fd, where we can take 
subunit-size d = I without losing any generality. At 
the stall force / = the filament delivers the maxi¬ 
mum work The free energy input to the 

filament in order to do this work is provided by poly¬ 
merisation and can be written as i^poiy = per 

subunit addition. Note that D monomers do not poly¬ 
merize, and hence there is no contribution in Apoiy due 
to D monomers. Finally we define the following quantity 
for a single filament 


a = Apoiy - W^diy = Mu/wt) - (5) 

The quantity a signifies how different the maximum work 
produced per filament is, as compared to the free energy 
input. Hence, a is analogous to the thermodynamic ef¬ 
ficiency of the system [^. In Figl2Kb) we plot the de¬ 
viation and the efficiency a for two microtubules, 
against the dissociation rate wd of D monomers. Inter¬ 
estingly, and a are correlated in the numerical sign. 
This gives us a hint that the violation of force additivity 
has something to do with the imbalance between the work 
produced and the free-energy input (i.e., the departure 
from equilibrium). To appreciate this interconnection be¬ 
tween A^^) and a, we proceed to investigate the particle 
fluxes of the T and D monomers. 

We calculate the particle fluxes when the A-filament 
system is at stall. In the simulations we separately 
keep track of the numbers of T and D monomers binding 
(or unbinding) at a filament-tip in a A-filament system. 
The fluxes are then calculated over a time window. The 
flux for T monomers per filament (Jr) is defined as the 
net change of T monomer numbers at any one filament- 
tip, divided by the size of the time window. Similarly, 
we also calculate the flux per filament for D monomers 
(Jd). In Fig. [21(c), we show these fluxes at stall. Al¬ 
though we have Jt + Jd = 0 (which is expected at stall), 
individually the fluxes per filament (both Jt and Jd) are 
non-zero signifying the non-equilibrium nature of the dy¬ 
namics. An important point to note is that, the fluxes per 
filament, at stall, decrease with the number of filaments, 
A, and tend to saturate. From this observation, we are 
tempted to make a hypothesis that a non-equilibrium sys¬ 
tem with a large number of filaments is closer to “equi¬ 
librium” in comparison to a single-filament system{see 
Appendix B for a crude entropy production argument 
for this hypothesis). 

With this hypothesis in hand, we now attempt to ex¬ 
plain why the numerical signs of A^^^ and a are corre¬ 
lated (Fig. EKb)). We first consider the case a > 0, when 
a single filament performs less work than the free-energy 
provided by the polymerization (see Eq. [5]). In this case. 
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FIG. 3. (a) Schematic diagram of a generalized random hydrolysis model with two-way switching (both T—^D, and D—^T). 

Different processes (shown in arrows) are discussed in the text, (b) Schematic depiction of a connected loop in the configuration 
space of a single hlament, within the model, (c) Deviation versus ut (the red curve), and a versus ut (the blue curve), 
for the generalized random hydrolysis model. The parameters are: wt = 2s~^,kTD = 0.3s~^, knr = 0.4s“^,u_d = 3s“^, and 
wd = ls“^. 


some energy is dissipated by the filament due to the in¬ 
ternal T—i-D transitions. As per our hypothesis, if the 
two-filament system is closer to equilibrium as compared 
to a single filament, the two-filament system should ex¬ 
tract more work by increasing the stall force per filament. 
To check this, in Figlljd), we show the stall force per fil¬ 
ament (i.e., the maximum work extracted per filament), 

fs^'^ /N as a function of the filament number N. The stall 
force per filament indeed increases with N for a > 0, and 
saturates near the net free-energy input (ln(u/i/; 7 ’)). In 
other words, as the number of filaments increases, the 
collective stall force per filament gets closer to the “equi¬ 
librium” value of la{u/wT)- This increase in stall force 
per filament makes /i^V^ > and in return gives pos¬ 
itive Thus, the a > 0 case correlates with > 0. 
Similar arguments can be given for a < 0 (see FiglJKd)), 
where a single filament performs more work than the en¬ 
ergy provided by polymerization. To bring the system 
closer to equilibrium, the two-filament system decreases 
the stall force per filament (/i^^/iV < /i^^). Hence, A^^^ 
is negative if a < 0. 

An interesting point to note in Figl2l(b) is that both 
A^^^ and a are zero exactly at wt = wd- This shows 
that T—>D switching (hydrolysis) is necessary to produce 
the phenomenon of non-additivity of stall forces. It is to 
be noted that the hydrolysis is always a non-equilibrium 
process as it is unidirectional. However, the condition 
Wt = Wd effectively corresponds to absence of switch¬ 
ing, since dynamically there remains no distinction be¬ 
tween T and D subunits. The filaments cannot “sense” 
their distinct presence as far as the force generation is 
concerned. Yet, the condition wt = wd does not imply 
a true equilibrium until we set the T^D switching rate 
to zero. Another way to possibly achieve equilibrium 
at stall is to incorporate the reverse switching (D—i>T) 
and allow the polymerization of D subunits. Although 
these additions are biologically unrealistic, we neverthe¬ 
less study such a model in the next section to explore the 
relevance of non-equilibrium dynamics for non-additivity 


of the stall forces. 


B. A generalized random hydrolysis model for 
filaments 

We make the random hydrolysis model (discussed in 
section IHI Ap more general and symmetric by allowing 
(i) D T conversion, and (ii) addition of both D and T 
monomers to the filaments. In this model (see Figl3](a)) 
both T and D subunits can bind to a filament with 
constant rates ut and ud^ respectively. When the fil¬ 
aments come in contact with the wall (see FigI3](a)), the 
polymerization rates decrease to ut (/) = UTe~^ and 
uoif) = UDe~^ in the presence of the force / (using the 
load distribution factor d = 1 for simplicity). The de¬ 
polymerization occurs with a rate wt if the tip-monomer 
is T, or Wd if it is D. Any randomly chosen subunit in¬ 
side a filament can convert either from T to D (with rate 
kro), or from D to T (with rate kox)- 

Within the general version of the random hydroly¬ 
sis model, we now proceed to show that the two-way 
switching (T —^ D, and D T) in general produces non¬ 
equilibrium dynamics that is embodied in the violation 
of the condition of detailed balance for the kinetic rates. 
For a single filament, as shown in Fig. [3Kb), we consider a 
loop of dynamically connected configurations. The prod¬ 
uct of clockwise and anticlockwise rates are urkToWo 
andunkDTWTj respectively. For the condition of detailed 
balance, that is, equilibrium, to be reached at steady 
state, the two products must be equal according to the 
Kolmo gor ov’s criterion (or the Wegschieder con¬ 

dition [ 93 ), which leads to 

Ut kTD Wd ^ 

-T-= 1- (6) 

Wt kDT ud 

If we fix the parameters wt = ‘2s~^,kTD = 

0.3s~^,kDT = 0 As~^,ud = 3s“^, and wd = then 

we would have ut = 8s“^ from the above equilibrium 
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condition (Eq. [Q. Though this criterion is written in 
terms of force-free rates, it is clear that using the modi¬ 
fied rates in the presence of the resisting force would not 
change the condition in Eq. El We now plot the deviation 
^(2) ^ 

versus ut in EiglSjc) ~ see the red 
curve (data from stochastic simulation). The plot quite 
interestingly shows that = 0 only at ut = 8s“^; 

otherwise it is nonzero. This clearly indicates that the 
phenomenon of non-additivity of stall forces is tied to the 
departure of the system from equilibrium. This can be 
compared with the arguments given in sectionlHl where it 
is shown that the filament models involving no switching 
exhibit equilibrium at stall, and consequently the rela¬ 
tionship fs^'^ = Nfs^'^ holds without any restriction. 

We can further relate the non-equilibrium effect of non¬ 
additivity of the force with the imbalance between the 
free energy input and maximum work output (as dis¬ 
cussed in section BlI Al) . Because, in the present model, 
the filaments can grow by adding D or T monomers, 
the free energy input to the system should depend on 
polymerization energies corresponding to both T and D 
monomers. Consequently, the partition function for a 
single filament is -I- (using ksT = 1), where ct 
and ejj are polymerization energies provided by T and 
D monomers, respectively. Hence, the free energy input 
to the filament is Epoiy = ln[e'^'^ -I- = h\[{uT/'WT) + 

{ud/wd)]] and the maximum work done by the filament 
is simply (using subunit size d = \). Following pre¬ 
vious section UlI Al we again define an efficiency-like pa¬ 
rameter for the current model as below 


a = T^poiy - = ln[(MT/ wt) + (ud/wd)] - (7) 


We see in FigEKc) that both a and are correlated 
in numerical sign, and they are non-zero every where ex¬ 
cept for a single point where the system is in equilibrium 
at stall (according to Eq. [6]). With this understanding 
of the connection between the non-equilibrium dynam¬ 
ics and non-additivity of stall forces, we study another 
model for filaments 31| in the next section, and show 
that the same ideas can be carried forward. 


C. A two-state model for filaments 

In the literature [Sllll , there exists a model which 
incorporates the detailed process of hydrolysis in a more 
coarse-grained way. In this model (Figl4l(a)), each fil¬ 
ament can switch between two chemical states I and 2, 
with switching rates ki 2 (from 1 to 2), and ^21 (from 2 to 
1). In the states I and 2, the filament has distinct depoly¬ 
merization rates wi and W 2 , respectively, and polymer¬ 
ization rates of ui and U 2 , respectively (see Figl3Ka)).If a 
filament bears the load (i.e., touches the wall), its poly¬ 
merization rate is modified to Ui{f) = Uie~^ {i = 1,2). 
For simplicity, we assume that the depolymerization rates 
are force independent. 


For the above two-state model, we first explicitly show 
that at stall the dynamics is non-equilibrium. As shown 
in FigUKb), we consider a loop of connected configura¬ 
tions for a single filament characterized by its length and 
state. In this case, Kolmogorov’s criterion reduces to 


Uiki2W2k21 = ki2U2k2lWi 
Ui _ U2 
Wi W2 


( 8 ) 


Following the procedure involvin g m icroscopic master 
equations, as described in Ref. [31|, we analytically 
obtain the single-filament and two-filament stall forces 


{fs^\ /i^^) In FiglUc), we plot A^^^ against U 2 (red 
curve) with fixed wi = 0.1s“^, ui = ls“^, and W 2 = 
0.8s“^. We find that A*^^^ ^ 0 for all U 2 , except for 
U 2 = 8s“^ (the value corresponding to the equality in 
Eq. [5]). This shows that the non-additivity of stall forces 
is tied to the departure from equilibrium. 

It is to be noted that in the two-state model, \n{ui/wi) 
is the polymerization free energy in state i = (1,2). 
Hence, if Pi is the probability of finding a filament in 
a state f, then at any instant the amount of free en¬ 
ergy that is transferred from the bath of monomers to 
the filament by addition of monomers of type 1 or 2 is 
T^poiy = -Pi ln(Mi/wi)-|-P 2 \n{u 2 /w 2 )- On the other hand, 
the maximum amount of work done by a filament against 
the applied force is = fs^'^ (monomer size being 

d = 1). Therefore, as defined in the previous sections, we 
can again define an efficiency parameter as: 


a = Fpoiy-W™®'’' = [Pi \n{ui/wi) + P2ln{u2/w2)]-fP, 

(9) 

where. Pi = fc 2 i/(fci 2 -I- fe 2 i) and P 2 = fci 2 /(fci 2 +^ 21 ). 
These probabilities follow from the fact that the de¬ 
tailed balance relation Pifci 2 = P 2 k 2 i holds at the steady 
state for a single filament, as intuitively evident from the 
FigUDa) (also see Ref M). along with the normalization 
condition Pi -|- P 2 = 1. In FigUKc), we plot a versus U 2 
(blue curve). We see that A^^) is closely coupled to a in 
numerical sign, and both are nonzero except at the point 
where the equilibrium condition (Eq. [8|) is satisfied. 

The effect of non-additivity of stall forces is not spe¬ 
cific to cytoskeletal filaments. Even the system of multi¬ 
ple molecular motors show such a behavior [s^, [Zlj . In 
the ensuing sections, we explore the connection between 
the non-additivity of the force with the non-equilibrium 
dynamics in the system of motors. 


D. Model of interacting motors by Campas et al. 

A model of multiple interacting motors pushing against 
a load has been recently proposed by Campas et al. [s^ 
. In this model (Fig. EJi-), motors walk along a one¬ 
dimensional lattice (lattice spacing d = \) and move by a 
single step forward (rate u) or backward (rate w). There 
is hard core interaction between the motors. The lead¬ 
ing motor alone bears the load, and hence its hopping 









FIG. 4. (a) Schematic diagram of single-filament 2-state model with switching between states 1 (blue) and 2 (red). Various 
processes are shown by arrows and corresponding rates, as discussed in the text, (b) Schematic depiction of a connected loop in 
the configuration space of single-filament 2-state model. The configurations are denoted by ordered pairs, whose first element 
is the instantaneous length and second element is the chemical state (1 or 2). (c)A^^^ (red) and a (blue) as a function of U 2 
for the 2-state model. Parameters are fci 2 = 0.5 s“^, /c 2 i = 0.5 wi — 0.1 ui = 1 'W 2 = 0.8 s“^. 


rates are modified to u{f) = and w{f) = 

where d is the force distribution factor. The hopping 
rates also change due to nearest-neighbour interactions 
- if a motor is adjacent to another one, its forward and 
backward hopping rates become u and w, respectively 
(see Fig. [IKa)). 


(a) 


u a/{f) u(f) 

,0, , ,Or-f, 


^ 0(f) 


(b) 


, 0 , 0 , 


i-1 i i + 1 i + 2 




a. 


i-1 i i-i-1 i-i-2 


•u/ 


U 


zi> 

IT 


,o, ,o. 


zv u 

,0, , ,Q, 

i-1 i i + 1 i + 2 


(C) 



u (s"'') 



FIG. 5. (a) Schematic diagram of the motor model proposed 
by Campas et al. [s^], where multiple interacting motors 
push against a cargo with a constant force / acting against 
their motion. Various processes are shown by arrows and 
discussed in the text, (b) Schematic depiction of a closed 
loop of dynamically connected configurations for the model 
shown in (a). (c)Deviation versus u within the model of 
Campas et al (data obtained from the exact formula given in 
[33|1. Parameters are specified in the text. We took the force 
distribution factor <5 = l.(d)Motors walking on a free-energy 
landscape [3^ . The effect of nearest-neighbour interactions is 
shown schematically. 


From analytical calculations and numerical simula¬ 
tions, the authors have found that the stall forces are 
not necessarily additive. We show here that the non¬ 
additivity is a manifestation of the non-equilibrium na¬ 
ture of the dynamics. To show this we apply the Kol¬ 
mogorov criterion, by making a closed loop of connected 
configurations as shown in FiglHKb). By equating the 


clockwise and counter-clockwise products of rates along 
the loop, we have 


u.u.w.w = w.u.u.w 
u u 


( 10 ) 


Exact analytical expressions of the single-motor stall 
force = ln(tt/rt)), and the two-motor stall force 

/o') 

f) ( = \t).\(uu/ww) + [u/w] — {u/w)] are derived in Ref 
[321 . If we put the equality of Eq. [TO] in the expression 
of ff‘\ we clearly see that the relationship = 2/i^^ 
follows. We further plot (Fig. EKc)) — 2/i^^ 

versus u, with fixed u = 20s“^,w = 5s“^, and w = ls“^. 
In this case, u = corresponds to equilibrium dy¬ 

namics (Eq. (To]), and we see that except for u = 4s“^, 
^ 0 everywhere. This corroborates our hypothesis 
that the additivity of stall forces is closely linked with 
the underlying equilibrium of the system, a point not 
recognized in Ref. [s^- 

We can also obtain the condition of equilibrium 
(Eq. IlOp from a simple thermodynamic argument by con¬ 
sidering the hopping processes of the motors on a free- 
energy landscape (see Fig. [SKd)). Due to the nearest- 
neighbour interactions, the shape of the free energy land¬ 
scape gets altered when the motors are adjacent to each 
other. An attractive interaction deepens the energy wells 
by an amount eint (dotted blue curve in Fig. [SKd)), and 
makes it hard for the particles to leave the position. This 
reduces the forward and backward hopping rates. On 
the other hand, the repulsive interaction makes the po¬ 
tential wells shallower (dotted red curve in Fig. EKd)), 
which makes it easy for the particles to hop forward or 
backward. In this case, hopping rates increase from the 
original values. When the motors are adjacent to each 
other, following Fig. [51(d), the hopping rates are given by 


w 

u 


= e 


e+fiint 


( 11 ) 

( 12 ) 


For a motor which is away from the other one, eint = 0, 
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and we further have 


u 

w 


= e 


(13) 


By rearranging the equations [TTl [12] and HS] we get back 
the same condition, as in Eq. (TO] which is a 

reflection of the equilibrium dynamics at stall. 


E. A motor model with multiple step-sizes 


It was found MM recently that dynein motors on a 
microtubule can take multiple step sizes - predominantly 
24-nm and 32-nm steps. This inspired us to cast a new 
model of motors walking with two distinct step sizes (see 
Fig. EKa)). A motor at a lattice site i can hop to any of 
the sites i + 1 (with rate iti), i + 2 (with rate U 2 ), i — 1 
(with rate wi), or i — 2 (with rate ^ 2 ). The leading mo¬ 
tor alone bears the applied force / and its forward rates 
are modified to ui(/) = Mie“^ and M 2 (/) = W 2 e“^^, while 
the backward rates are assumed to be force-independent. 
Unlike the model in the previous section (Section 1111DD 
there is no explicit attractive or repulsive interaction be¬ 
tween the motors. 

Proceeding in a similar way as described in section 
IIII B[ we first derive the condition for equilibrium, by 
considering a closed loop of connected configurations as 
shown in Fig. EJb). Following Kolmogorov criterion we 
get 


u\.W2 = U2-wl 



(14) 


If we take ui = 80s“^,r(;i = 8s“^, and W 2 = ls~^, then 
U 2 = 100s“^ is the value corresponding to the equilib¬ 
rium condition (see Eq. fUl) . In Fig. EKc), ver¬ 

sus U 2 is plotted from our simulation results (the red 
curve), where we see that indeed = 0 only for 

U 2 = 100s“^. Thus, one can again associate the force 
inequality ^ N to the violation of the detailed 
balance condition. 

Just like the models of filaments discussed in the pre¬ 
vious sections (Sections IIII A[ IIIIBI and IIII Cl) , we can 
show that the effect of force-additivity for motors is re¬ 
lated to the underlying energetic imbalance. The motor 
model with multiple step-sizes looks conceptually similar 
to the generalized random hydrolysis model (see Section 
IIII Bl) . At any instant, a filament can grow by adding 
T or D monomers in the generalized random hydrolysis 
model; whereas in the current model a motor can move 
forward (or backward) by taking steps of sizes d or 2d. 
These steps of size d and 2d involve different amounts of 
work done by a single motor at stall (/i^^d and fs^^2d re¬ 
spectively). By contrast, addition of a D or T monomer 
to a single filament leads to the extraction of the same 
amount of work (/i^^d). Hence, for the definition of the 
efficiency parameter, a, we choose to take into account 


the energy imbalance that results from a single step of 
size d. For a single motor taking a step of unit lat¬ 
tice size, we can write the free-energy supplied by ATP 
molecules as F = ln(Mi/r(;i) and the maximum work done 
as (considering the lattice spacing d = I). 

Thus, we can define an efhciency-like quantity as 

a = F-W'^^^ = ln{ui/wi)-fi^l (15) 

We do not claim that this is a unique definition of a for 
the system. One may certainly come up with some other 
definition of a, for example, 

a= ^(In—-2/^(1)), (16) 

2 W2 

to quantify the excess/deficit of energy supply to the sys¬ 
tem. Note that both the definitions of a (Eas.lTblandflBll 
essentially capture the energy imbalance per unit step. 

We plot a (from Eq.[Tn|) and A^^) versus U 2 in FigEIJc), 
and see that both a and A^^^ are indeed correlated in 
numerical sign. Moreover, both of them are zero only 
when the equilibrium condition (Eq. [T^ is satisfied. The 
same finding can be derived using the other definition 
of a, that is, using Eq. [T6| instead of Eq. [T5| (data not 
shown). This further strengthens our point that non¬ 
additivity of stall forces is a manifestation of underlying 
non-equilibrium stall dynamics. 


IV. DISCUSSION AND CONCLUSION 

Collective force generation by filaments or motors 
has been theoretically studied by many researchers us¬ 
ing various m odel s in specific contexts [H, m (IB- 
Ill) [13, EB, [73 llOOl . Il0l| . However, a broad picture ex¬ 
plaining the cooperative effects in stall (maximum) force 
generation is still missing. In this paper we have pro¬ 
vided a theoretical framework to understand and predict 
the cooperative effects in the maximum force generation 
by multiple motors or filaments, for a broad class of mod¬ 
els. It is now appreciated, at least theoretically, that the 
stall force of individual cytoskeletal filaments or molec¬ 
ular motors, when they push together against some ob¬ 
stacle, is not additive in general [U, EB, [53|. In this 
paper, we have provided several pointers to show that 
this non-additivity of the stall forces (/i^^ ^ Nfs^'^) is a 
manifestation of the underlying non-equilibrium nature 
of the dynamics at stall. 

Our study suggests that if the Kolmogrov criterion for 
kinetic rates is satisfied then one does not need a detailed 
non-equilibrium calculation to obtain the stall force of 
multiple filaments/motors. The same result can be ob¬ 
tained from a simple equilibrium statistical mechanics 
calculation. Moreover, even if the Kolmogrov criterion 
for kinetic rates is not satisfied for a system, our effi¬ 
ciency parameter a qualitatively predicts that the coop¬ 
erative effects in the stall force for multiple filaments (or 
motors) is either enhanced or decreased as one increases 
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FIG. 6. Multiple step size motor model: (a) Schematic diagram of the motors moving with two distinct step sizes. The kinetic 
processes (shown in arrows) are explained in the text, (b) A closed loop of connected conhgurations for the model with one motor, 
(c) Deviation versus U 2 (the red curve), and a versus U 2 (the blue curve). The parameters are ur = 80s“^,U)r = 
and W 2 ~ 



FIG. 7. Schematic of the (a) continuous and (b) dis¬ 
crete two-state Brownian ratchet model [102| . 1 1031] . (a) A 

spontaneous motion is expected when the ratio of transition 
rates is far from the equilibrium value given by the de¬ 

tailed balance condition, (b) Discrete version of the two-state 
ratchet model with effective transition rates a{x), uja{x), 
ujb{x) and tj s( 2 :).(c)Energy landscape corresponding to dis¬ 
crete two-state model. 


the number of filaments (or motors). For a class of mod¬ 
els discussed in this paper, we clearly see a correlation 
between the numerical signs of a and the force deviation 
We would like to point out that the Kolmogorov or 
Wegschieder criterion for detailed balance has been used 
earlier in the context of multiple growing filaments [2^ 
to describe the linear scaling of the stall force with the 
number of filaments. We, however, used the Kolmogorov 
criterion to systematically probe a range of models and 
derive the conditions that must be obeyed by their re¬ 
spective kinetic rates, if the systems are expected to be 
in “equilibrium” at stall. More importantly, the concep¬ 
tual advantage gained by couching this problem in the 
language of thermodynamic equilibrium is that, if any 
similar system is not in equilibrium at stall, then it is 
not at all guaranteed that the stall forces will be additive 
with the system size (number); in fact, non-additivity is 
the more likely outcome. 


To illustrate the above point further, we take a con¬ 
crete example of a two-state motor model from the ex¬ 
isting literature [H, [H, [tI, [t^ (Fig. [T]), which de- 
mostrated the stall force non-additivity ^ Nfs^\ 
Our analysis indicates the same conclusion without the 
need for any detailed simulation or theoretical calcula¬ 
tions. We identify that the violation of detailed balance 
in the transition rates between the states is a require¬ 
ment for the spontaneous motion in this two-state Brow¬ 
nian ratchet model (see Fig. EKa)). As a result, even 
for a single motor at stall, although the mean velocity 
flux of the motor is zero, the individual velocity fluxes 
in states-I and state-2 (with potentials Ui and U 2 , re¬ 
spectively) for the motor are not independently zero. In 
fact, only by observing the energy landscape one can say 
that, at stall, the mean particle flux in state-I and state- 
2 would be positive and negative, respectively. This is 
because, the mean particle flux in state-2 is zero in the 
absence of any resisting force (flat potential, see FiglTKa)) 
and would become negative with the application of a re¬ 
sisting force. To make the mean flux of the overall sys¬ 
tem zero, the effective mean flux in state-1 has to be 
positive under stall condition. This non-equilibrium at 
stall for one motor manifests itself in the form of non¬ 
additivity of stall forces in the presence of multiple mo¬ 
tors even though the only interaction between them is 
self-exclusion (FiglTKa)) [^[7l,[^. In terestingl y, a dis¬ 
crete version of this two-state model |in2l . llOdj . shows 
additivity of stall forces for multiple motors with only 
self-inclusion interactions (see Fig. EKb)). This results 
from the fact that this discrete model for a single mo¬ 
tor can easily be mapped to a biased random walk with 
only one track [s^ - the motors effectively move on a 

tiled energy landscape, e = In (see Fig. (Tl), which, 

OJA B 

as argued in Section II, indicates that the motor can be 
interpreted to be in equilibrium at stall. 

We can generalise the observations noted above and 
propose that the stall behaviour of passively interact¬ 
ing, processive motors greatly depends on the topology 
of the mechanochemical migration path of the motors. 
Consequently, we generally propose that motor mod¬ 
els analogous to the ones pioneered by Kolomeisky and 
Fisher [J with a single track for motor migration 
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will show additivity of stall force for multiple, passively 
interacting motors. On the other hand, the Brownian 
models with multiple tracks demonstrably exhibit 
non-additivity of stall forces for multiple, even passively 
interacting motors [zl, [t^. To the best of our knowl¬ 
edge, the general conceptual framework that compares 
and contrasts the behaviour of these two major classes 
of molecular motor models, has not been provided before 
and is also one of the main contributions of this paper. 

To summarize, collective force generation in biofila¬ 
ments and molecular motors typically involve multi-step, 
complex internal processes and a variety of interactions 
between the individual entities as well as the source of the 
resisting force. In this paper we, however, have demon¬ 
strated that, even with very simple internal dynamics 
and also in the absence of any attractive or repulsive in¬ 
teractions between individual components, if a system 
of molecular motors or filaments is not in equilibrium at 
stall we can expect non-additivity in their collective force 
generation. To establish this result with reasonable cer¬ 
tainty, we have analyzed, multiple seemingly disparate 
models, which nevertheless exhibit a common theme of 
non-equilibrium dynamics at stall leading to this coop- 
erativity. The formalism developed in this paper should 
provide a general thermodynamics based framework with 
which to perform primary interpretation of experimental 
and theoretical results relating to collective force genera¬ 
tion in biofilaments and molecular motors, before exam¬ 
ining the system related specifics. 
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APPENDIX A 

SIMULATION METHOD 

We have simulated our models using kinetic Monte- 
Carlo algo rithm als o known as the Gillespie algo¬ 
rithm, flM flo^ , for the calculation of stall force 
of analytically unsolvable models . Now we elaborate 
the exact method we have used to simulate “general¬ 
ized random hydrolysis model’ for n = 3 microtubule 
filaments. We start time evolution of the model sys¬ 
tem from a fixed initial length of Iq = 2000 monomer 
at t = 0 for all the microtubules. To obtain the time 
at which the next event will occur, we sum all possible 
event rates, oq, and generate a random number p be¬ 
tween 0 and I. Now r = — In ( ^ ) is the time of 

ao \p J 

the next event. To determine which event k will occur 


at time r, we generate another random number p 2 be¬ 
tween 0 and 1 and find k which satisfies the condition 
Si=i o,i < ®oP 2 < Si=i Now we repeat this pro¬ 
cesses till our system reaches a steady state. For n = 3 
microtubule, we have taken time to reach steady state as 
t = 1000s“^. Now at this point we reset our time to zero 
and start monitoring the position of largest microtubule 
filament with time. We plot this time evolution data and 
from the slope of the plot we get the velocity of the tip of 
largest filament. We perform this simulation for various 
forces and obtain corresponding velocity of the barrier. 
We then draw the force velocity curve. The point at 
which this force velocity curve cuts the velocity axis is 
the stall force. To further verify the stall force, we cal¬ 
culate particle flux of the monomers for that force. For 
the calculation of flux, we just mark one filament and 
observe the number of monomers being added to that 
filament (-1-1) and removed from that hlament(—1). Now 
we divide this count of monomers by the time window 
for which this bookkeeping has been done. If this flux is 
close to zero(> 10“^) then we take that point as the stall 
force. 


APPENDIX B 

CONNECTION BETWEEN SUBUNIT FLUX 
AND ENTROPY PRODUCTION FOR MULTIPLE 
BIOFILAMENTS AT STALL DESCRIBED USING 
THE RANDOM HYDROLYSIS MODEL OF 
SECTION [iHAl 

In the following text, we provide a simple, and perhaps 
crude, connection between the subunit flux jr^N) per 
filament, presented in Fig. [2j: for the random hydrolysis 
model of Sec. lIII Al and the corresponding coarse grained 
entropy production s{N) per filament for N Hlament at 
stall. 

Entropy production per hlament is given by [tI, Il02j| : 

s = vf + rAfj,, (17) 

where /, v, Ap and r are the external force, mean ve¬ 
locity, chemical potential for the subunit exchange and 
subunit exchange hux, respectively. At stall condition, 
u = 0, and hence, the entropy production is induced 
only by subunit exchange hux. For the simple T-D ran¬ 
dom hydrolysis model described in Sec. IIII Al the effec¬ 
tive subunit huxes (per hlament) are jx and jo, respec¬ 
tively. Moreover, at stall, since the total material in- 
hux, on an average, into the system of hlaments is zero, 
jT = —jo = j- Since, in this model, the only subunit 
on rate is U for the T subunits, the chemical potentials 
corresponding to T and D inhux can simplistically be 
written as 

A 1 U 

A/tt = In-, and 

Wt 

A 1 U 

Apo = In-. 

wo 


( 18 ) 
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Using Eqs. [17] and (TH the entropy production at stall 
condition can now be re-written as 


= jT{N)A^T + jD{N)AflD 


= j{N)\n^. 

Wt 


(19) 


Thus, a very simple argument proposes that entropy pro¬ 


duction per filament at stall condition is proportionate to 
the subunit flux per filament in the system of multiple 
hlaments. 
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